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Abstract. We present 850 /im imaging polarimetry of the W51A massive star forming region performed with 
SCUBA on the JCMT. From the polarimetry we infer the column- averaged magnetic field direction, projected 
onto the plane of the sky. We find that the magnetic field geometry in the region is complicated. We compare the 
field geometry with 6 cm and CS J=7-6 emission and determine that the magnetic field must be relatively weak 
and plays a passive role, allowing itself to be shaped by pressure forces and dynamics in the ionised and neutral 
gases. Comparisons are drawn between our data and 1.3 mm BIMA interferometric polarimetry data, from which 
we conclude that the magnetic field must increase in importance as we move to smaller scales and closer to sites 
of active star formation. 
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1. Introduction 

The W51A cloud consists of a complex of Hn regions 
situated approximately 7-8 kpc away. It is one of the 
most active and luminous sites of massive star formation 
in the Galaxy with a ll the associ ated signatures: masers 
(Lep panen et al |l998 ; Zh ang fc Ho 1995; Genzel & Downes 



1977j ) , Hn regions (Scott l978|), bright & dusty cores (Jaffe 



et al. 1984[ Genzel et al. 1982) and n earby supernova rem- 
nants (Koo & Moon 1997a , 1997b ). Radio maps show a 
number of Hn regio ns, th e most luminous being G49.5-0.4 
(see e.g. Mehringer 1994 ) which itself is dominated by the 
sources W51d and W51e. The former is directly associated 
with its i nfrar ed counterpart IRS2 (Goldader & Wynn- 
Williams 1994 ) and while W51e and the infrared source 
IRS1 are also closely associated, this section of the cloud is 
more complex. A cluster of compact Hn regions are found 
to the east of W51e/IRS1, whose radio emission is dom- 
inated by the sources W51e2 and el. The submillimetre 
emission is more closely associated with this cluster rather 
than with IRS1 and the peak of emission correlates with 



the position of W51e2 (Jaffe et al. 1984), suggesting that 
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W51e is a more evolved source while W51el and e2 are 
only now undergoing star formation. This is corroborated 
by the detecti on of infalling gas associated with W51e2 by 
Zhang & Ho (|1997| ). 

Here, we present 850 /im polarimetry of the W51A re- 
gion revealing the magnetic field structure within it. There 
are a number of methods which allow for mapping the 
magnetic field geometry, and each has its own associated 
problems. By far the most promising method is through 
polarimetry of emission from dust grains which have been 
aligned (it is necessarily assumed) by interaction with the 
local magnetic field. Previous polarimetric observations of 
W51A show very low polarisation (Kane et al 1993) consis- 
tent with a non-detection. However, these measurements 
were made with large beams (~ 30") within which a sig- 
nificant degree of depolarisation can occur. For instance, 
Dotson et al. ( [2000 ) present far-infrared polarimetry for 
W51A at 100 /im taken with the KAO. Their results show 
low polarisation at the cores and higher degrees of polar- 
isation in regions with more extended emission. Higher 
resolution and more sensitive measurements are clearly 
necessary to detect significant polarisation and to deter- 
mine any s tructure in the pola risati on pattern (e.g. see Lai 
et al. 2001 and Momose et al. 2001 ). In their mid- infrared 
spectropolarimetric atlas of young stellar objects, Smith et 
al. ( 200C ) measure 6% and 3% polarisation towards W51d 
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in absorptive and emissive components, respectively, mea- 
surements made within 5" beams. 

The importance of magnetic fields in star formation 
has become increasingly more apparent in recent years as 
computer models have become more sophisticated (e.g. 
Ostriker et al. 2001). Magnetic fields are believed to play 
an important regulatory role during cloud collapse and ac- 
cretion onto the protostar, as well as providing the driving 
and collimation mechanism for outflows. However, in or- 
der to comprehend the magnetic field's role properly, an 
empirical understanding of its behaviour is necessary. 

In this paper, we present 850 /im polarimetry not only 
of the dense cores in W51A but also of the surrounding 
more diffuse emission between these cores. In this fash- 
ion we investigate the relation between the magnetic field 
at the sites of star formation and the cold, expansive en- 
velopes that surround them. By comparing this informa- 
tion with other studies we may also hope to understand 
what shapes the magnetic field. For instance, are the kine- 
matics of gas and dust, as governed by gravitational col- 
lapse and the pressures within Hn regions, dominant, or 
is the energy density in the magnetic field of sufficient 
magnitude to control and regulate the dynamics? 

2. Observations and Reduction 

The data for W51A were observed on the night 6 October 
2000 on the James Clerk Maxwell Telescope (JCMT) sit- 
uated close to the summit of Mauna Kea, Hawaii. The ob- 
servations were made with the SCUBA camera (Holland 
et al. 1999) together with a polarimetry module in the 
beam mounted on the entrance window of SCUBA. The 
polarimetry optics consist of a rotating quartz half-wave 
retarder ahead of a fixed wire-grid analyser (see Greaves et 
al. 2000 for a full description of the polarimeter). To pro- 
duce a fully sampled map at 850 /im, 16 offset positions 
of the array are needed on the sky (Holland et al.). The 
resultant image is referred to as a "jiggle map" . For these 
observations of W51A an east-west chop throw of 140" at 
7.8 Hz was used. Complete jiggle maps were taken at spe- 
cific positions of the half-wave retarder, each separated by 
22.5°. This therefore requires 16 jiggle maps to complete 
one cycle of the retarder. Six such cycles were observed 
towards W51A. 

The SUR F (SC UBA User Reduction Facility; Jenness 
& Lightfoot 1998|) data reduction package was used to 
flatfield, rebin and sky noise correct the SCUBA data. 
Excessively noisy bolometers were identified from noise 
measurements conducted on the observing night and 
"switched-off" so that they take no further part in the 
data reduction. Sky subtraction with SURF is used for 
eliminating sky noise in SCUBA data (Jenness et al. 1998[ ). 
Sky noise can be introduced, for example, by variations in 
the sky transparency between the separate positions of 
the secondary, which make up the jiggle pattern, in a sin- 
gle exposure. To correct for this noise bolometers devoid 
of any significant flux which can be attributed to point 
sources, or to extended emission, are chosen. The signal 



from each jiggle position in these bolometers is then used 
to calculate and correct for any temporal variation in the 
sky signal. To preserve the total flux in the map, the mean 
background signal which is removed in this correction is 
then added back. Inspection of the final image showed 
that the weakest signal is an order of magnitude fainter 
than that of the diffuse emission. This weak signal being 
~ 2a above the sky noise, suggests that conditions were 
stable enough over the timescale of the chop cycle that sky 
background was efficiently removed by the chop. The very 
small, residual sky emission would have to be quite highly 
polarised to have any effect on our results and conclusions. 
With the sky noise removed, the data are then rebinncd 
onto a 3.09" pixel grid, representing half the spacing of 
the jiggle pattern, to form the final images. 

Before rebinning the images and calculating the po- 
larisation, the main lobe instrumental polarisation (IP) is 
removed from each bolometer using a look-up table of el- 
evation dependent IP values. This software also corrects 
for the parallactic angle at the time of observations by al- 
tering the value of the effective retarder position on the 
sky. Following Greaves et al. (2002) and Matthews et al. 
(2001), contamination due to sideband lobe instrumental 
polarisation is estimated to be negligible (i.e. within er- 
rors) for our data. 

The /, Q and U Stokes parameters were computed for 
each cycle of the retarder using the POLPACK polarime- 
try data reduction software (Berry & Gledhill 1999 ). For 
each cycle, images taken at equivalent positions of the 
retarder are averaged together (e.g. positions at 0°, 90°, 
180° and 270° equivalently measure the same angle on the 
sky) . The result is a stack of four images where each pixel 
in the stack samples the modulated signal at angles close 
to 0°, 22.5°, 45° and 67.5° to which the following curve is 
fitted to extract the Stokes parameters (Sparks & Axon 



1999): 



I' = -(7 + e[Qcos20 + [/sin2</>]) 



(1) 



where I' is the expected intensity, t is the throughput of 
the wire-grid analyser, e is the polarising efficiency of the 
analyser, and <f> is the effective retarder position angle after 
correction for the parallactic angle. 

After combining the data from the six cycles, the 
Stokes parameters are binned into 3x3 pixel bins (equiv- 
alent to 9.3"). Finally, the polarisation is debiassed and 
thresholded so that only polarisation vectors associated 
with positive flux (/ > 0) and small polarisation error 
(SP < 0.5%) are included in the final data set. The bulk 
of our polarisation vectors lie in the 1-3% range which, 
when the maximum polarisation error of SP — 0.5% is 
taken into account, translates to an uncertainty for the 
position angles of A6 ~ 28.6°/ct p - 15 - 5°. 

Finally, flux calibration was performed by multiply- 
ing the data by a "flux calibration factor" (FCF) of 414 
Jy/beam/Volt. This FCF has been shown to be repeatably 
accurate to ~ 10% for all SCUBA calibrations (Jenness et 
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Table 1. Brightest Hll regions and their positions in 
the W51A cloud. Positions obtained from the SIMBAD 
database. 



Object 



RA (J2000) Dec(J2000) 



W51d (IRS2) 
W51e (IRS1) 
W51el 
W51e2 



19 23 39.9 
19 23 41.9 
19 23 43.8 
19 23 43.9 



+14 31 08 
+14 30 36 
+14 30 25.9 
+14 30 34.8 



al. 2002). The value used here has been modified to ac- 
count for the loss of throughput due to the polarimeter 
(Greaves et al. 2002). 

3. Results and Discussion 

3.1. Submillimetre polarimetry 

The final results are tabulated in Table 2 and shown in Fig. 
|l|. The greyscale shows the submillimetre emission which 
is dominated by two bright sources. To the east is the core 
which was mapped by Jaffe et al. (1984) at 400 /im and 
is associate d with the W51el and e2 compact Hll regions 
(Mehringer 1994 ). The submillimetre core to the north- 
west is associated with the W51d Hll region an d wit h 
the IR source IRS2 (Goldader & Wynn- Williams |l994|) . 
Both co res a re known sites of maser activity (Genzel & 
Downes 1977 ). Fainter and more extended emission is also 
seen, occupying the region between the two cores, and an 
extended arm of emission is seen to reach out to the south- 
west from IRS2. 

Fig. |] shows polarisation vectors overlaid onto the sub- 
millimetre emission. These vectors have been rotated by 
90° to depict the column-averaged direction of the mag- 
netic field projected onto the plane of the sky. This field 
direction towards the W51el/e2 core is at a position an- 
gle of ~ 100-110°, consistent with recent interferometer 



polarimetry at 1.3 mm with BIMA (Lai et al. 2001 ). In 
an atlas of mid-infrared spectropolarimetry of young stel- 
lar objects, Smith et al. fl2000| ) observed W51d(IRS2) and 
deconvolved emissive and absorptive components to the 
spectropolarimetry. The emissive (E-vector) component 
was at 36° ± 3° implying a field direction in the hotter 
regions at or near the Hll region (there is a clear [Nell] 
line in the spectrum) where the 10 jxm emission peaks 
of 126°. The absorptive component of the E-vector, and 
therefore equivalent to the B-vector, is at a position angle 
of 136°, a difference of ~ 40°when compared to our mea- 
surement of ~ 0°. Being absorptive, the dust responsible 
for the mid-infrared polarisation is also probably largely 
responsible for the submillimetre emission. This difference 
in the inferred magnetic field directions is most likely due 
to structure in the magnetic field being sampled to differ- 
ent optical depths at the two wavelengths (and the differ- 
ent beam sizes to some extent). 

Scatter plots of the polarisation data (Fig. ||) show 
that the degree of polarisation has a maximum of ~ 5% 
and also that the polarisation signal-to-noise, up, ranges 
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Fig. 2. Scatter and histogram plots of the polarimetry pre- 
sented in Fig. |l[ The top panel shows the degree of polar- 
isation variation with the logarithm of intensity. The bot- 
tom panel shows the variation of op, with polarisation. No 
points exist below the dotted line due to the data-clipping 
of points with polarisation errors SP > 0.5% (see Sect. 2). 

up to a maximum value ~ 20. The top panel of Fig. || 
shows that across the region the polarisation and inten- 
sity are generally anti-correlated, a trend seen by others 
(e.g. Matthews et al. |200l| ; Davis et al. feOOCp . A number 
of physical mechanisms have been attributed to explain 
this effect. A decrease in grain alignment efficiency as a 
function of increasing depth into a cloud is often cited 
as the cause (e.g. Lazarian et al. 1997). This is because 
most known mechanisms of grain alignment are believed 
to decrease in efficiency at the physical conditions com- 
mon to cold, dark clouds. On the other hand, small scale 
field structure induced by field lines being dragged and 
twisted under the influence of gravitational collapse, for 
instance, are capable of producing the depolarising effect 
if the telescope beam is larger than the scal e leng ths over 
which these effects occur (e.g. Aitken et al. 1997 ). Zhang 
et al. ( 1998j ) present evidence which suggest that cores 
within IRS1 are indeed collapsing and rotating. However, 
as a caveat to this it should be noted that this trend of low 
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Fig. 1. SCUBA polarimetry of the W51A region. The intensity image is taken from our data and is plotted on a linear 
scale. The source to the east of the field is the submillimctrc core associated with W51el and e2, and the source to the 
west is W51d. The polarisation vectors have been rotated by 90° and plotted such that their position angles delineate 
the direction of the magnetic field averaged along the column and projected onto the plane of the sky. The white 
dots represent the positions of the brightest Hn regions in the field (see Table [l]). They are, from right-to- left: W51d 
(IRS2), W51e (IRS1), W51el and W51e2 (el is south of e2). Right ascension and Declination are at epoch 2000. 



degrees of polarisation for high submillimctre intensities 
is not universal as can be seen in recent submillimetre po- 
larim etric imaging of the NGC7538 region (Momose et al. 
2001 ) where the degree of polarisation remains relatively 
high towards the cores. 

It is also worth noting that polarisation nulls can be 
seen to the north-east and south-west of W51el/e2, at 
positions lying almost perpendicular to the field direction 
defined at the core. These null vectors lie at symmetric 
positions, that is, they are both equispaced and on the 
same line through the source centre. The data have been 
checked and there is no apparant systematic reason for 
these vectors to be lowly polarised apart from the fact 
that the polarimetry at these positions has low signal-to- 
noise. A physical reason may be that at these positions the 
field becomes wound up such that the polarisation drops to 
very low values (P < 0.5%) and the polarised intensity is 
averaged out within the relatively large JCMT beam. We 
choose to briefly discuss these features here because simu- 



lations do show that such polarisation patterns are possi- 
ble forpinched and twisted field morphologies (Aitkcnct 
al. $Q0l\) as well as for helical fields (Fiege & Pudritz p000| ) , 
although the latter deals with linear structures and is not 
really applicable to discrete sources. Furthermore, such 
'submillimetre polarisation nulls' either side of a bright 
core have been previously observed towards NGC7538 
IRS1 (Momose et al. 2001). In both of these cases, the 



field pattern appears very much like the familiar polari- 
sation disks seen in scattered light, at near- infrared and 
optical wavelengths, around sources with optically thick 
disks (e.g. Lucas & Roche 1998). Higher spatial resolu- 



tion observations would clearly help us to understand the 
physical basis behind these polarisation nulls, as we note 
that there are nulls present in our data at other positions 
in the region. 

An immediately striking feature of the polarimetry as 
a whole is that the magnetic field threading the region 
does not appear uniform in direction or degree of polarisa- 
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tion. It is interesting to note that the polarisation vectors 
which are positioned central and to the south in Fig. ^ 
lie closely aligned to the projections of the Galactic plane 
and the large scale Galactic magnetic field onto the plane 
of the sky (PA - 20°; Matthewson & Ford |197l| ). In their 
study of water masers in the core of W51e2, Leppanen et 



al. (1998) found the masers' proper motions also aligned 
in this direction. Muson & Liszt (1979) find a large scale 
stream of molecular gas at this same postion angle pass- 
ing through the centre of the W51A region. Their CO 
data shows a complex of velocities with the bulk motion 
at 58 km s -1 , the streamer velocity at 70 k m s -1 an d cold 
absorbing gas at ~ 66 km s _1 . Koo & Moon ( 1997b ) deter- 
mined that 13 CO emission is associated with the streamer. 

Clearly, this position angle of ~ 20° bears some signif- 
icance for this region. Those vectors which closely follow 
this position angle, may be tracing the large-scale Galactic 
magnetic field which follows the Galactic plane and may 
also be dragged and stretched by material in the streamer. 
The field then becomes disrupted when it encounters the 
dense, star-forming cores. Towards W51el/e2 and W51d 
the field appears perpendicular to this global direction, 
suggesting that some degree of field dragging and winding- 
up must have occured during the initial collapse of these 
cores. 

The arm of submillimetre emission which extends to 
the south-west of W51d also has the field associated with 
it perpendicular to the Galactic plane and the arm itself. It 
appears that the collapse to form this elongated structure 
has occured along the magnetic field lines, but this does 
not seem probable if those field lines are connected to field 
lines to the north and to the south of the structure, where 
the field apparantly joins the Galactic Plane. It could be 
that the gravitational collapse has dragged the field lines 
into this pattern. 

3.2. Comparison with 6 cm and CS J=7-6 emission 

We have obtained a 6 cm emission map of the W51A 
region, downloaded from the NCSA Astronomical Data 
Image Library. The data were originally acquired at the 
VLA and published by Mehringer (1994). Also, raw data 
from a CS J=7-6 raster map of the W51A region were re- 
trieved from the JCMT archive (operated by the CADC) 
and reduced. Fig. ||shows contour plots of the 6 cm and CS 
emissions overlaid onto our submillimetre measurements. 

The 6 cm map traces free-free emission from ionised 
gas and the eastern peak, W51e, coincides with the IRS1 
source seen in the infrared. The extended structure of the 
6 cm emission seems to correlate with the magnetic field 
direction quite well. The extended emission as a whole 
gives the impression of an envelope around the western 
edge of the submillimetre core and, from the polarisation 
pattern, the magnetic field appears to be influenced by 
the presence of this ionised gas. There is an extension to 
the north-east of W51e which follows the magnetic field 
lines and may form part of the molecular streamer seen by 



Mufson & Liszt (197E). In all, the 6 cm emission possesses 
a degree of curvature which may be a consequence of the 
expansion of the Hn region being impeded by the dense 
submillimetre core. If one follows the inferred magnetic 
field direction from the north-east through to the south- 
west, then one finds that around the location of W51e the 
field is indeed curved. 

Being a density tracer, the integrated intensity map of 
the CS J=7-6 emission looks very similar to the submil- 
limetre emission. However, the contour map presented in 
Fig. H shows the integrated flux between 59 and 63 km s _1 , 
and at these velocities, a bridge of material between the 
two cores becomes apparant. Channel maps of the CS 
data show that there exists a small velocity gradient of ~ 
3 km s _1 across the bridge. This suggests a flow along the 
bridge and between the two cores, although this evidence 
may be better interpreted in terms of a rotating structure. 
Proper dynamical modelling, beyond the aims and scope 
of this paper, is necessary for a better understanding of 
this structure. 

The relation of the magnetic field with the material 
responsible for the CS J=7-6 emission is not as clear as 
for the 6 cm emisssion. The most intriguing feature in the 
CS emission is clearly the bridge structure between the 
two cores, a feature not pronounced in the 850 /im emis- 
sion. From the W51el/e2 core, the bridge extends to the 
north-west. Here the magnetic field seems to lie perpen- 
dicular to the bridge, although it could easily be more in- 
fluenced by the molecular streamer in the Galactic Plane 
(see Sect. 3.1) or the nearby W51e Hn region. Further 
along, the bridge then curves around to be more latitudi- 
nal; the magnetic field fans out and is now parallel with 
the bridge. If there really is material flow along this struc- 
ture, then it may be that the magnetic field is being pulled 
and stretched along the bridge at this position. 



3.3. An upper limit to the magnetic field strength 

As well as the morphology of the magnetic field, it is also 
important to try and obtain estimates of the magnetic 
field strength. Methods that have been used in the litera- 
ture have included the Zeeman effect (e.g. Crutcher 1999 ), 
which gives an indication of the strength of the field com- 
ponent along the line of sight, and also the method of 
Chandrasekhar & Fermi (1953) which gives the strength 
of the field component along the plane of the sky through 
measurement of the dispersion of polarisation position an- 
gles. The latter method is applicable to our data set but 
would not return sensible results due to our large beam 
and also because the large dispersion of position angles we 
find across our field rend ers this method futile (see Sect. 
6 of Ostriker et al. 2001 ). However, we are able to place 
quite reliable upper limits to the magnetic field strength 
through consideration of the Virial Theorem. For a self- 
gravitating molecular cloud with negligible external pres- 
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Fig. 3. Submillimetre polarimetry image of Fig. [I] overlaid with contours of 6cm emission (left panel; Mehringer 1994 , 
courtesy of ADIL) and CS J=7-6 emission (right panel). The presence of the W51e Hn region associated with IRS1 
seems to be influencing the magnetic field. The CS emission (integrated between 59 and 63 km s _1 ) reveals a bridge 
of dense gas between W51el/e2 and W51d (IRS2). Once again, the magnetic field appears to be influenced by this 
structure, especially the northern portion of the bridge. Right ascension and Declination are at epoch 2000. 



sure the virial equilibrium of that cloud can be written as 
(McKee et al. |1993|; Glenn et al. |1999Q: 



-3GM 2 



5R 



= 2 x 0.27MAv 2 + 0.1B 2 R 3 



(2) 



For W51e2 we use a mass of 80 (Zhang et al. 
1998j ), Av ~ 1 km s" 1 (Young et al. |1998[ ), and a radius 



10 m, measured from the 



for the cloud of R 
1.3 mm map of Lai et al. (2001), to obtain a value of 



B up ~ 3.5 mG. For W51d we use a mass of 400 Mq 
(Zhang et al. 1998] ), Av ~ 1 km s^ 1 again, and a radius 
for the cloud of R ~ 5 x 10 15 m, measured from the 
2 mm map of Zhang et al. (1998), to obtain a value of 



0.6 mG. 



Since we know that these clouds are in the process of 
collapse and forming stars, these values form useful up- 
per limits to the value of the magnetic field strength. The 
value for W51e2 is consistent with values obtained us- 
ing the method of Chandrasekhar & Fermi (~ 1 mG; Lai 
et al.). Crutcher (199E) presented a statistical analysis of 
Zeeman measurements from a number of molecular clouds 
(but not W51) and found line-of-sight field strengths of the 
order of 0.5 mG for clouds with densities of 10 5-6 cm -3 . 



4. The Magnetic Field in W51A: Just How 
Important Is It? 

In Sect. 3 above, we describe how the magnetic field fol- 
lows the curvature of the extended 6 cm emission and is 



also influenced by the bridge of CS J=7-6 emission be- 
tween the two submillimetre cores. The apparant curva- 
ture of the field and the manner in which the field aligns 
with the velocity gradient in the CS bridge can be inter- 
preted as evidence that the magnetic field is being shaped 
by the pressure and dynamics of the gas. Similar magnetic 
field behaviour has been seen in recent SCUBA polarime- 



try of NGC 7538 (Momose et al. [200 1[ ) where the blue- 
and red-shifted outflow lobes from IRS 11 are observed to 
have quite a high degree of curvature. The magnetic field 
inferred from the polarimetry follows the same curvature. 
However, one would expect that if the field were of suf- 
ficient strength, magnetic tension in the field lines would 
resist such action and return them back into an ordered 
and linear configuration. 

Towards the core of the W51el/e2 region we find that 
the magnetic field direction is oriented roughly perpen- 
dicular to the long axis of the extended envelope. This is 
clearly suggestive that the material collapsed along the 
field lines and then fragmented to form the cores now 
present. This interpretation is advocated by Lai et al. 



(2001) who found a very uniform field direction across 
the two cores (with the benefit of 3" resolution from their 
intcrfcrometry data) and interpreted this to mean that the 
magnetic field plays a dominant role. However, rotational 
motions have been detected t oward s the collapsing cores 
of W51el and e2 (Zhang et al. 1998 ), and their rotational 
axes are between 40-90° from that of the magnetic field. 
This misalignment is evidence that the magnetic field has 
had no influence in determining the rotational axes for 
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these clumps, evidence that their fragmentation and col- 
lapse from the common envelope has been guided accord- 
ing to physical parameters other than the magnetic field. 

At the small scales studied by Lai et al., their reason- 
ing may well hold true but the evidence presented here 
suggests that the magnetic field in W51A does not play 
a dominant role in determining where gas flows. Firstly, 
across the 2' field-of-view we find a large dispersion of 
polarisation position angles (see Fig. [|). Recent MHD 
modelling of a turbulent molecular cloud (Ostriker et al. 
2001) has shown that, as expected, clouds with weak fields 
(/3 ~ 1, where (3 is the ratio of the squares of the sound and 
Alfven speeds and parameterises the relative strength of 
the magnetic field; see Eq. [|) produce polarisation maps 
with vectors which are less organised and ordered than 
those with stronger fields ((3 ~ 0.01). This is consis- 
tent with what we are observing in W51A (see Fig. 23 
in Ostriker et al.). 

We may analyse this further by investigating the di- 
mensionless parameter, (3. This parameter is defined as 
(Ostriker et al. [2001]): 



/3 



(3) 



where c s is the sound speed, va is the Alfven velocity, Bo 
is the magnetic field strength and p is the mass density. 
Using a sound speed of of = 3fcT/m, we can rewrite the (3 
parameter as the ratio of the gas to magnetic pressures: 



nkT 



(4) 



where n is the particle number density, T is the temper- 
ature and k is Boltzmann's constant. Taking estimates of 
density and temperature from the literature we find that 
n ~ 4 x 10 5 - 10 6 cm" 3 and T ~ 20 - 35 K (see e.g. 



Lai et al. 2001 and Jaffe et al. 1984). We can now write 
the (3 parameter specifically for the W51A region as: 



(3 ~ 10- 8 /B 2 



(5) 



In their simulations, Ostriker et al. defined a 'strong 
field' condition when (3 = 0.01 which produced well cor- 
related and ordered polarisation maps. This translates to 
a field strength of Bo ~ 1 mG in W51A. Conversely, 
their weak-field case has (3 — \ which translates to a 
field strength of B ~ 0.1 mG. 

We therefore find a consistency between our data, 
which infers a weak field on large scales, and that of Lai et 
al., who infer a strong field on smaller scales. The magnetic 
field has an important role to play in the W51A region 
as one approaches the dense cores where star formation 
is active. However, as one moves to the more expansive 
envelope and to regions between star forming cores, the 
magnetic field appears weaker and plays a more passive 
role, allowing itself to be shaped by the pressure forces 
and dynamics of the gas. 



5. Conclusions 

We present 850 pm polarimetry of the massive star form- 
ing region in W51A. The presented data has relatively 
high polarimctric signal-to-noise, which allows us to infer 
a magnetic field structure which is not uniform across our 
2' field-of-view. 

We have compared our data with a VLA 6 cm emis- 
sion map as well as with a CS J=7-6 map obtained from 
the JCMT archive. The comparison shows a correlation 
between the magnetic field geometry and the ionised and 
molecular gas. This indicates that the magnetic field is 
shaped by the gas dynamics and not the reverse, and 
therefore has a relatively weak magnetic field strength (pa- 
rameterised by (3 ~ 1). 

However, Lai et al. (2001) have found from their 



1.3 mm interferometer polarimetry that at smaller scales 
and towards the star forming cores the field lines are reg- 
ular and ordered, implying that the magnetic field is com- 
paratively strong (f3 ~ 0.01) and plays a significant role 
in determining the gas dynamics - although evidence sug- 
gests that further fragmentation and collapse probably oc- 
curred independent of the magnetic field. Nevertheless, 
when one looks at larger scales and further comparisons 
are made between the global field morphology and the gas 
dynamics, this picture breaks down. 

Our data has shown how important it is to attempt 
to correlate magnetic field structures with the physical 
conditions of the (ionised and neutral) gas. In this way, one 
can develop a better understanding of the role magnetic 
fields play in star formation. 
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